Quantum dots with Rashba spin-orbit coupling 
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We present results on the effects of spin-orbit coupling on the electronic structure of few-electron 
interacting quantum dots. The ground-state properties as a function of the number of electrons 
in the dot N are calculated by means of spin density functional theory. We find a suppression of 
Hund's rule due to the competition of the Rashba effect and exchange interaction. Introducing an 
in-plane Zeeman field leads to a paramagnetic behavior of the dot in a closed shell configuration, 
and to spin texture in space. 
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Spin-related phenomena have attracted great attention 
recently as they are the key ingredient in the emerging 
field of spintronics M . Among these, spin-orbit (SO) cou- 
pling mechanisms in semiconductors provide a basis for 
device applications, and a source of interesting physics, 
especially in systems with reduced dimensionality. Trans- 
port through chaotic quantum dots in the presence of SO 
interaction has been studied both experimentally O and 
theoretically H ||, while the effect of SO coupling on the 
spin lifetime has been investigated in Ref. ||. Here, we 
are interested in how the electronic properties of few- 
electron quantum dots, such as the addition energy]^, |t]] 
or the spin properties of the dot ground state M, are af- 
fected by Rashba SO interaction || M|. These questions 
are interesting from the theoretical point of view for the 
following reasons. First, the Rashba effect has a different 
form than the usual SO coupling term in real atoms. Sec- 
ond, the tunability of the Rashba effect JTT|, p"2[ |l3| allows 
dramatic SO effects to occur in quantum dots with few 
electrons; in real atoms this requires heavy nuclei, and 
hence a more complicated electronic structure. 

We start by describing the physics of a quasi zero- 
dimensional system with Rashba SO at the non- 
interacting electron level, providing analytical results for 
the single-particle spectrum when the SO coupling can 
be treated as a perturbation. Then, we introduce the 
electron-electron interaction in the framework of Spin 



Density Functional Theory (SDFT)Q. This allows us to 
study how the addition spectrum is modified by varying 
the strength of SO coupling. Studying the spin properties 
of the many-particle ground state, we find a suppression 
of Hund's rule, when the SO coupling can still be treated 
as a perturbation for the single-particle problem. For 
higher strengths it affects the single-particle spectrum so 
strongly that it gives rise to a completely different addi- 
tion spectrum. The introduction of an in-plane magnetic 
field, leads to a paramagnetic behavior of the dot in a 
closed shell configuration, and to spin texture in space. 

Quantum dots are often realized by lateral confinement 
of a two dimensional electron gas (2DEG) obtained in a 
heterostructure. Due to the lack of inversion symmetry 
along the growth direction z of the heterostructure 0, nfl, 



the electrons in the 2DEG are subject to the Rashba spin- 
orbit coupling Hamiltonian 



(cr x Py - CTy Px) 



(1) 



The strength of the SO coupling, here denoted as k so , 
can be tuned by changing the asymmetry of the quantum 
well via externally a ppl ied voltages, as shown in several 
experimental studies |Tl], [1^, [ll| . 

It is interesting to study the effect of the SO coupling 
term Eq. (Q) on the quantum mechanics of a quasi zero- 
dimensional system [{l5|. To this end, we consider a two- 
dimensional quantum dot defined by a parabolic confin- 
ing potential 



V can f(x,y) = — w (x + y ) 



(2) 



Thus, the single-particle Hamiltonian in the effective- 
mass approximation reads 



H= Px+Py +Vco ^ xy ^ +Hso 

2m 

In the absence of SO coupling the eigen-energies are 



(3) 



E$ =huj (Af+1) 



(4) 



with M being a non-negative integer. A degenerate sub- 
space <Sm of dimension Dm = 2(M+ 1), where the factor 
2 is due to spin, is associated to each energy E^ . 

We will now treat -ff so as a perturbation. This is valid 
as long as k so l u <C 1, where l u is the oscillator length 
\Jhj (mw). We obtain for the second-order eigen-energies 



E 



E<$ + M^oU 2 [2(t - I) - (M + 1)] , (5) 



where i = 1, • • • , M + 1, and a = ±1 is the quantum 
number relative to cr z , i.e. the spin projection along z. As 
the single-particle levels will play an important role in the 
following, we show an example of the low-energy part of 
the spectrum calculated using Eq. (||) together with the 
results of numerical diagonalization in the upper panel 
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of Fig. [|. From the perturbative treatment the following 
conclusions can be drawn: 1) Each degenerate level E$ 
is split in M sublevels, each of which is double degenerate 
due to Kramers theorem; 2) Spin rotational invariance 
is broken but still the eigenstates are (to this order in 
perturbation theory) eigenstates of a z . From conclusion 
1) we can infer that SO coupling changes the addition 
spectrum of the dot, while conclusion 2) tells us that the 
Rashba effect will not influence the lifetime of the eigen- 
states of a z . For values of k so for which perturbation 
theory breaks down, the eigen-energies are of course still 
grouped in Kramers-degenerate sublevels [as shown in the 
inset of Fig. 0(a)], although it can happen that different 
sublevels have almost the same energy. 

To introduce the Coulomb interaction between the 
electrons we use spin density functional theory (SDFT), 
in the local density approximation . We write the 
Kohn-Sham equation []16| in spinor notation 



2m 



conf 



coul 



Ha, 



(6) 



where V^ ou i is the Hartree potential, V exc h- C orr the 
exchange-correlation potential, 1 the identity matrix in 
spin space, and is a two-components spinor. The spin- 
dependent exchange-correlation potential is introduced 
following Ref. H]. ft is in general non diagonal in spin 
space (it becomes diagonal in the basis that diagonalizes 
the spin-density matrix). For the exchange-correlation 
energy we use the Tanatar and Ceperly parameterized 
formp7|; in particular for the case of partial spin polar- 
ization we use the interpolation scheme of Refs. |8|, p"4[ . 
We solve Eq. (||) self-consistently by discretizing it in real 
space. The addition energy (also called capacitive energy 
by some authors^) is defined as E^dd = (i(N+l)—[i(N), 
where n(N) is the chemical potential for the dot with N 
electrons, i.e. the energy needed to add the N-th. elec- 
tron to the system containing already N — 1 electrons. 
We compute the chemical potential by means of Slater's 
rule 1 18 , in order to minimize numerical errors due to 



double differentiation. 

Now, we focus on realistic dots obtained in an InAs 
heterostructure, where the Rashba effect can be quite 
large Q. We use for the electron effective mass the value 
m = 0.022 too, with mo being the free-electron mass; and 
for the dielectric constant e = 14.6 eo, being eo the one 
of vacuum. 

In a quantum dot without SO interaction we expect 
peaks in the addition energy when the number of elec- 
trons N equals a magic numbers, i.e at those integer val- 
ues which correspond to a closed shell configuration. For 
a parabolic dot the first magic numbers are 2,6,12, see 
Eq. (0) and the discussion below it. Besides these peaks 
for N coinciding with a magic number, some additional 
peaks are expected for a number of electrons correspond- 




FIG. 1: (a) Low-energy part of the single-particle spectrum 
for the dot, calculated both by perturbation theory (dots) 
and by numerical diagonalization (squares) for ksolu = 0.2633 
[this value corresponds to the data for fc so = 0.01 nm _1 for the 
addition energies of the realistic dot shown in panel (b)]. The 
label p is an index that enumerates the eigenstates in order 
of ascending energy. Inset: Low-energy part of the single- 
particle spectrum calculated numerically for k S olu> = 0.7896 
[this value correspond to fc so = 0.03 nm _1 in panel (b)] In this 
case the SO coupling dominates the single-particle spectrum, 
(b) Addition energy vs number of electrons in the dot for 
different values of the SO coupling strength. In this figure 
and in the following ones the dot is defined by a confining 
potential of strength fiu) — 5 meV. 



ing to a half-filled shell due to Hund's rule. In this sit- 
uation the electrons in the half-filled shell have parallel 
spins to gain exchange energy. For the parabolic dot un- 
der consideration the first Hund's-rule peaks are located 
at N = 4, and N = 9. Both these kind of peaks can 
be seen in Fig. 0(b), for the case when no SO coupling 
is present (filled dots). Switching on the SO interaction 
leads to a change in the single-electron levels of the dot, 
and such a change is reflected in the addition energy. This 
effect can be seen in Fig. 0(b), where the addition energy 
is plotted for several values of the SO coupling strength, 
ranging from a situation where the perturbative treat- 
ment Eq. (Q) is still valid, to one where the SO coupling 
dominates the single-particle spectrum. Due to the pres- 
ence of Kramers-degenerate sublevels in the free-electron 
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spectrum, peaks tend to be present for even number of 
electrons in the dot. The fact that a degenerate level .Efi 
is split in M Kramers-degenerate sublevels (in the pertur- 
bative regime), leads to a suppression of the Hund's rule: 
in a half-hlled level to maximize the total spin the elec- 
trons should be allocated one per sublevel, but this has an 
energy cost equal to the sublevel splitting; if the sublevel 
splitting is larger than the gain in exchange energy, then 
the spin polarization for a half-filled level is suppressed. 
This is indeed what we see by analyzing the spin prop- 
erties of the SDFT ground-state wave-function. At this 
point it is important to stress that due to the tunability 
of the SO coupling strength O, |l2| , [ll|, it is possible to 
investigate experimentally the effect of the Rashba term 
on the addition spectrum of few-electrons quantum dots 
(addition spectra were measured by Tarucha et aZ.0), 
and the transition from weak (k so l M <C 1) to strong SO 
coupling. 

We now investigate the effect of an in-plane magnetic 
field. Due to the fact that the system is invariant under 
rotation around the z-axis, we can choose the direction 
of the in plane magnetic field arbitrarily without losing 
any generality. We introduce a magnetic field B along x, 
which does not affect directly the orbital motion, but cou- 
ples to the x-component of the total spin, giving rise to a 
Zeeman term, H z — huj z S x /2, where S x = 53;=i n a x\ 
and 7ko z = /iB5* B, with /ie being Bohr's magneton, and 
g* the g-factor. We consider now a dot in a closed-shell 
configuration, namely we take N = 2 and N = 6. In the 
case of vanishing SO coupling and in the independent- 
electron picture, such a system does not respond to the 
in-plane magnetic field for Zeeman splitting smaller than 
the level splitting (u> z < uj). The situation changes when 
the Rashba term is introduced; the ground state of the 
dot exhibits now some net-spin polarization. In the upper 
panel of Fig. ||, the average value of 5a, is plotted vs mag- 
netic field, showing how the system gets magnetized even 
in a closed-shell configuration due to the interplay of SO 
coupling and Zeeman splitting. The average values of 5 
and S z remain equal to zero. In the lower panel of Fig 
the variation of the ground state energy with magnetic 
field is plotted vs magnetic field. It shows a decrease with 
increasing field (which is well fitted by a parabola), yield- 
ing a positive susceptibility \ = —d 2 E/dB 2 . Thus, we 
can conclude that the dot in a closed shell configuration 
exhibits a paramagnetic behavior. This is in contrast to 
what happens in real atoms, where a closed shell gives a 
diamagnetic response due to orbital degrees of freedom 
(Larmor diamagnetism) [[l9[ , while in our case the Lar- 
mor term is suppressed by the two-dimensionality of the 
dot. This paramagnetic behavior [^0| is due to the single- 
particle eigen-states (see below), but it persists when the 
electron-electron interaction is present (it is enhanced by 
it), as shown in Fig. ^. 

In the limit of ksol^ui / lj z -C 1, and w z < w, it is pos- 
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FIG. 2: (a) Average values of S x vs the strength of the 
in-plane magnetic field, computed by means of SDFT, for 
two closed shell configurations (N — 2 and N — 6) and for 
different strength of the spin-orbit interaction. Without the 
Rashba term, (S x ) would be zero. 

(b) Variation of the ground state energy AE = E(B) — E(0) 
vs the strength of the in-plane magnetic field, computed by 
means of SDFT. The ground state energy shows a quadratic 
dependence on the in-plane magnetic field: E(B) = E(0) — 
AB 2 , with A > 0. The susceptibility is positive and the 
system is paramagnetic. 



sible to obtain a perturbative expression in H so for the 
single-particle eigen-energies: 



Ema,<j 



E 



(o) 

M 



huj z 
- — 

2 

, ,2 



huj 



i 



" (^so^u;) 



l +tr (2i-i)_£ , (7) 

OJ 



where i = 1 • • • M + 1, a = ±1 is the quantum number 
relative to the projection of spin in the direction of the 
magnetic field, i.e a x , and E^ are the energies given in 
Eq. (^). In the independent-electron approximation, the 
energy of a closed shell E CS (M) = J^. EM,i,a reads 



E CS {M) 



2(M + l)E {0) 



M 



-(M + i)huj{k so L) 2 



l 



• (8) 



Expanding Eq. (@) in uj z /uj we get for the magnetic- 
field dependent part of the closed-shell energy — (M + 
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FIG. 3: Spin density projected along x (a), y (b), and z (c), 
computed by means of SDFT, for a dot with N = 6 elec- 
trons, fcso = 0.02 nm -1 , and Tuo z = 1 meV. In (a), for ease 
of visualization, we plot —n x instead of n x . Spatial integra- 
tion of n y and n z gives zero, yielding a zero average for the 
corresponding total-spin components. 



l)%uj{ksoluj) 2 (tjJ z /ui) 2 + 0[(uj z /lj) 4 }, which explains the 
parabolic behavior seen in Fig. ^. From Eq. (^) we get a 
paramagnetic contribution to the susceptibility due to a 
closed-shell 



XM 



(S*MB) 2 



2(M + lV 3 (fc s X) : 



(9) 



which in the limit of w z /oj <gc 1 is just a positive constant: 

Xa//(.9>b) 2 = 2(M + l)(fc s X) 2 /(M- 

It is interesting to have a closer look at the spin-density 



for the magnetized dot. In Fig. || the projections of the 
spin density along the x-, y-, and z-axis are shown for a 
dot containing six electrons in the presence both of SO 
coupling and of a Zeeman field. As it is clearly visible, 
the system shows spin texture in space, this is due to 
the fact that no common spin-quantization axis exists 
anymore (a similar situation occurs in quantum wires 
with strong spin-orbit coupling pl[.) 

In conclusion, we have studied the effect of Rashba 
spin-orbit interaction on the addition energy, and on the 
spin properties of a few-electron quantum dot by means 
of spin density functional theory. In particular, we have 
found a suppression of Hund's rule, for small fc so values, 
for which perturbation theory in H so still holds. An ad- 
ditional in-plane magnetic field (Zccman field) leads to a 
paramagnetic behavior of the dot in a closed shell con- 
figuration, and to spin texture in space. 
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